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Abstract 

Using an asymptotic solution of the A/-impurity thermodynamics of a 
dilute s-d system, the impurity energy and impurity heat capacity AC(T) 
are derived for dilute magnetic alloys with spin 1/2 and spin 3/2 im- 
purities. The parameters which enter AC are adjusted to fit exper- 
imental data on impurity heat capacity of CuCr and (Lai-xCe^) AI2. 
Agreement is satisfactory for CuCr, at temperatures below IK, and good 
for (Lai-zCe^) AI2. The magnitude of theoretical AC(T) agrees with ex- 
periment and does not require scaling as in previous s-d theories. Nonlin- 
ear dependence of AC(T) on impurity concentration has been accounted 
for the first time. 



1 Introduction 

Existing theories of anomalous thermal behaviour of dilute magnetic alloys 
(DMA) (e. g. Refs. PMI1I1GI1IS1I71IE1GI]) have provided only partial quanta- 
tive explanation of these anomalies. Most of these investigations were restricted 
to an s-d system consisting of a single magnetic impurity interacting with the 
electron gas. 

Significant progress was made by Andrei et al. [UGlj, who solved the s-d ther- 
modynamics for an s-d system with impurities treated as indistinguishable parti- 
cles. Their heat capacity and magnetization curves, after scaling by an adjusted 
factor, provide a good fit to experimental data on (Lai-zCe^) AI2. 

The disadvantage of almost all existing theories of DMA is linear depen- 
dence of resulting thermodynamic functions on impurity concentration c. The 
dependence of DMA experimental data on c is more complex, e.g. the temper- 
ature at which DMA resistivity minimum occurs is proportional to c 1 / 5 [T] and 
normalized DMA impurity heat capacity AC/c is not constant in c [10) . 
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A different approach to the thermodynamics of DMA, in which the termody- 
namic functions do not require rescaling and depend on c nonlinearly, was pro- 
posed by the author in Ref. [TT] . The starting point in this approach is the first 
quantization s-d Hamiltonian H^ n ' M ' describing n electrons and M distinguish- 
able, arbitrarily positioned impurities with single-component spins interacting 
via a locally smeared s-d interaction U(R a — fi)S za a Z i. 

(M n \ 
H o l) +.9 2 EE U &* ~ ^)S za o zl , (1) 
a=l i=l / 

where denotes the antisymmetrizer with respect to electron variables with 
indices i = 1, . . . , n, 

fc2 n 
i=l 

R a denotes the position of ath impurity, S za its spin operator and f*j, a Z i the 
corresponding quantities of the ith electron. The considerations of Ref. [TT] 
remain valid for any sufficiently regular function U depending on \R a — fi\. 
It was shown in Ref. [TTJ, that in the dilute limit of small c (d-lim) the s-d 

interaction is separable and the system's free energy per electron f(Hj n,M \f3) 
is equal asymptotically to that of n noninteracting electrons with the 1-electron 
Hamiltonian 

hi^H^v) =W' M) &V) + \M{i 2 V 2 % (3a) 

where 

M 

(£, „) = _ g ^ _ „) £ UP*™ (3b) 



and M impurities described by the Hamiltonian 

M ^ 

T 



M M 

C (0=flV^E 5 -+ V E ■ ( 4 ) 



Ua in Eq. ([3b]) is the multiplication operator by U(R a —ri) and rj(£) — £— /2(£), 
with 

h ® = -^ s *kz ' {B) » := TT Z e -^ > (5) 

whereas £ is the solution of the equation 

e = /i(/ 2 (0) + /2(0 (6) 

with 

F^B^ := A {n) (b {1) ® l {n -^ A {n \ 
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which minimizes the mean-field free energy per electron f{h^ M \f3), where 

h^ M \z, v ) = ht M \z, v ) + ftPJ(0. (8) 

Asymptotic equivalence of Hj 71 ' and h^ n ' M ^ is expressed by the equality 
lim d-lim/(iJ I ( "' M) ,/3) = lim &-\iuif{h {n ' M \ (3) = lim f{A^H^\[3). (9) 

n — >oc n — >oo n — >oo 

The question arises, to what extent is the mean-field s-d Hamiltonian h^ n,M ^ 
capable of providing a reliable theory of DMA thermal behaviour. Here a partial 
answer to this question is found by showing that impurity heat capacity curves 
of two CuCr alloys [TU] with c = 21.7 ppm, 51 ppm and (Lai-zCe^A^ [U] with 
x = 0.0064 can be satisfactorily explained in terms of /j("> M ) thermodynamics. 

To this end, the 1-particle equilibrium density operator p^ of a quantum 
gas, in a field of randomly positioned wells, is first analysed in Section [2] Such 
operator appears, e.g. in the equation for the fugacity z: 

Tr zp {1) fa + zpW) 1 =n (10) 

where 



p^ := exp 



It is shown that, in the low-temperature regime, such a gas behaves effectively 
like a system of free particles at a temperature higher than that of the real 
system. 

In Section[3]the mean-field impurity energy At/ S _d of the s-d system, relative 
to that of the free electron gas, is derived for impurity spins 1/2, 3/2. Here the 
crucial question is the form of f\. An approximate Sommerfeld-type expansion 
of /i(£) is found, which is subsequently applied in the simplest truncated form, 
viz., = b + with 6q,&i € M 1 . As a consequence, A[/ S _d depends 

on five parameters: bo, b±, g, M , AT, M denoting the number of impurities in 
the considered s-d subsystem of the molar s-d system and AT > the shift in 
temperature scale due to random interactions unaccounted for by the s-d Hamil- 
tonian |(T]|. 

The expression for the mean-field heat capacity 

AC(T) = 

was calculated analytically and the parameters bo, b\, g, M, AT adjusted to 
obtain directly the best possible fit of mean-field AC(T)/c with experimental 
data, without additional rescaling procedures. 

Variation of AC/c with c [TU] requires different values of b , &i, g, M, AT for 
different impurity concentrations. In this manner, for the first time, nonlinearity 
of AC(T)/c with respect to c has been accounted for. 
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1-particle density operator of a quantum gas 
in a field of randomly positioned wells 



The Hamiltonian 



(n,M) 
I 



gives only a simplified account of the interactions 
present in an s-d system. It does not contain terms representing the Coulomb in- 
teractions between electrons and the screened Coulomb pseudo-potential J2 a ^Cu 
of each electron in the field of impurities. In order to investigate the effect of 
such potential on the electron gas in a DMA, let us examine the 1-particle den- 
sity operator of a gas of negatively charged particles interacting with randomly 
positioned positively charged ions in a region A. 
Let 

* (1, = -^+£<#2 ( 12 ) 



denote the 1-particle Hamiltonian. Suppose Uc is a sufficiently regular func- 
tion of \R a — r\, so that the integral kernel of p^ — exp [— /3/iW] admits 
the Feynman-Kac representation [T3HH], viz., 



p(r,r') = / d/i- r -,(u>)exp 



p 



(13) 



where Uc a (r) — Uc(R a — f)- Due to randomness of ion positions R a , p(r,r') 
is equal to its space average (^p(r, r ')/ A over these positions [TUGS]. In order 



to evaluate \p{r, r')j , let us note that Uc a is negative- valued and has a finite 

minimum [17]. Suppose it admits a Taylor expansion in the neighbourhood of 
this minimum at fo a and that [3 is sufficiently large. Then 



p(r,r>) ) =|A| 



-M 



M 

J] J d 3 R aP (r,r>) 



(14) 



can be evaluated by expanding Uc{R a — t) up to second order and applying 
the method of steepest descent. Evaluation of the average (Ti"4"|) thus amounts 
to calculating the integral 



P 3 

Y[ J d 3 R a exp -(3u - ^u 2 J ^ (R aq 



roo 



UJq(s)) d.S + 



where u = U c (R a -foo^j, 



"2 



q=l 



u's(r 



(15) 



roc 



This is done in Appendix El 
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One obtains for sufficiently large (3 

[p (?,?)) A =Zi I d// ( ., ; lw) / exp[-0Mu o 

dAidA2dA3 



9=1' 



where 



-Mu 2 J KO) " \f ds 

/^XI(M-I) 



(16) 



Z l = |A| 



M 3/2 



The integral kernel (fT6|) thus represents, up to a constant factor, the kernel of the 
canonical density operator of a particle oscillating around a point A, integrated 
over all positions A £ ffi 3 . In other words, 

|A|- M J] f d 3 R aP ^ = Zi exp[-/3M Uo ] / d 3 Aexp [-/3^(A)j (17) 



where 



ft 2 



/lose (A) = —ij-A-, 

2m 2 



-Mu 2 (f-Xj 
The integral kernel of exp -/3ft££(A)] is [H HI EEB] 
1 



Posc(r,r') = £ 3 exp 
with 



a 

2 



(r-\j +(p-\j +bff-\\(p-Xj 
(tanhifi) 172 - 2 



(18) 



(19) 



2V7rsinh|n 



, a = or coth il , b 



sinh SI 



y 2 = %- 1 ^m^ , n- n2a2 



id 



The integral over R 3 on the rhs of Eq. (fl6|) is thus a product of three 1-dimensional 
integrals 



dA exp 



a — b 



1 

— a 
2 



exp 



(x - A) 2 + (x' — A) +6 (x - A) (a/ - A) 
7 (a + 6) (s - x'f 



(20) 



each of which is equal to the integral kernel of the operator [T31 Q3] 



,-b 



V27rf exp [-iolb] = 



2tt 

Va 2 - b 2 



exp [-*o7b] 



(21) 
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where 



One finds 



with 



t Q T =t(u 2 ,P) (~-£s) =t(u 2 ,(3)H ( ol ) 



t{u 2 ,0) = 6~ x tarih(8P) , 8= — a 2 

2m 



(22) 



Combining Eqs. (|T6|). (|20|). ([21]) , (|22|) . one obtains for large (3 the asymptotic 
equality 



p(r,r')) = Z Q (u Q ,u 2 ) (2irt a ) 3/2 exp 

A 



1 

2h 



(23) 



where 



Z (u ,u 2 ) = |A| 



|(M-1) 



Thus 



-M f 

= Z exp 

/ A 



M 3/2^3 



2tt 



exp [— /3Mu ] . 



-tH, 



(i) 



(24) 



In the low-temperature regime, {p^) A is thus equal, up to a normalization 
factor, to the canonical density operator of a free particle at an effective tem- 
perature (ke^r 1 > (ke^r 1 = T. 



3 Mean-field impurity heat capacity 

Formulae for s-d system's mean-field energy AC/ s -d and heat capacity AC, rela- 
tive to that of free electrons, will be now derived, using the Hamiltonian h^ n ' M \ 
and compared with experimental data on impurity heat capacity of CuCr [lOj 
and LaCeAl 2 [15]. 



3.1 CuCr alloys 

According to Monod et al. [TO] the spin oftheCr 3+ ions in CuCr alloys equals 3/2. 
The impurity expectation energy in these alloys therefore equals 

(25) 



<-Amp — y l imp / ( M ) 



-M</ 2 (0 + -Mg 2 + AMg — 



2 y y e -^9 2 cosh (3/3g£y/n) + cosh (0g£y/n) ' 



G 



/i, defined by Eq. (|7|) is a linear function in the simplest approximation (jl5|) : 
/i(0 = 6 +6i£ (Appendix|B]). Then according to Eqs. (O, (gt]), the interaction 
energy of electrons with the Hamiltonian h e n ' M (£,??) equals 



■hi"-^^) \ (26) 
= -M»/ 2 (0 (6o + 6i/ 2 (0 + •■■)■ 

^From Eqs. ©, (f2l>l) , (f2"6")l one obtains the following expression for AU s -d of 
a system of n electrons and M impurities with spin 3/2: 

AC/ s _ d - (7 imp - Mn.m) (b + 6i/2(0 + ■•■) + M</ 2 (0 ~ ^Mn/|(0 (27) 

where £ is the minimizing solution of Eq. ^ with /1 given by Eq. (I45D . 

The number of conduction electrons per host atom in CuCr equals one, and 
hence the number of these electrons per impurity rt\ — c _1 . 

The n-electron, Af-impurity s-d system will be now treated as a subsystem 
of a sample S containing one mole of impurities. Then n = Mri\. In terms 
of m, 7 = VMg, the energy AU m = 6.022 • 10 23 Af -1 A{7 s _ d of such sample 
with spin 3/2 impurities, divided by c and expressed in joules, equals 



c 



- x A[/ n 




-4/5 7 2 M 1 cosh (3^^) 



(28) 



e-WM" 1 cosh (3/?7CV^T) + c «sh (P^^/m) 
-ImV/KO - M 2 m/ 2 (0 (6 + 6i/ a (f) +•■•)) 602.2 ■ 160.2 M"V 
with 

= 7 te-WM- 1 sinh (3^y^T) + sinh (/3 7 £Vm) 
M4j M^nT e- 4 ^ 2 ^- 1 cosh (3/? 7 £V™T) + cosh (/^f^/nl) 

if 7, 60 are given in V eV. 
Eq. (|6|) takes the form 

/s(0=£ (29) 

with/ 3 = 6 + (6i + l)/ 2 . 

The excess heat capacity of 5, relative to that of pure Cu, equals 



(31) 



where 

dT dT \ d£ 

AC(T) depends on the solution £ of Eq. ([29]) and on the parameters 60, 61, J,M. 
An additional parameter AT > 0, equal to the shift of the temperature variable 
of AC, proves necessary. The origins of this shift are explained in Section O 
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Alloy 


X 


C 


71 1 


6i 


fe [VeV] 


T[VeV] 


M 


AT [K] 


CuCr 




51 
10" 


10° 
51 


-461 


i.09 ■ iezf. 
«i_ 


0.091 


248500 


1.05 


CuCr 




212 

To 7 


10 7 

212 


-631 


1.01 -i^ 


0.086 


36 • 10 4 


0.78 


(LaCe) Al 2 


64 


4 
1871 


3742 

3 


-13101 


nr^ 


31 


58 


0.39 



Table 1: 



The best fitting graphs of c _1 AC(T + AT) vary with M. Moreover, the 
graphs obtained for M 3> 1 provide much better agreement with experiment 
than those for M = 1,2. The graphs of c~ 1 AC(T + AT) for CuCr with c = 
51ppm (M = 24850) and c = 21.2 ppm (M = 360000), plotted in the same 
units as in Ref. [10], are depicted in Figs. [TJ [2j The corresponding values of the 
remaining parameters bo, b\, 7, AT are given in Table [T] 

The solution £(T) of Eq. (|29[) is unique in both cases and minimizes f{hS n ' M ',fi). 
The graphs of £(T) are plotted in Fig. H 

It appears that higher order terms of the expansions (|4"2"|) , (|45[) should be 
included in order to improve the quality of mean-field AC/c. 

It is worth noting that since the best fitting values Mf of M are much smaller 
than A — 6.022 • 10 23 , the sample S can be viewed as consisting of magnetic 
domains, each containing Mf impurities with a definite favoured impurity-spin 
orientation, which differs, in general, from one domain to another. Existence of 
such domains in some magnetic materials has been established experimentally 
(e.g. Ref. [20]). 



3.2 (Lax_.Ce,) Al 2 

Experimental data on AC/c of (Lai-zCe^) AI2 alloys are presented in Ref. [12] . 
According to Refs. 12, 21j, a typical Kondo effect, without any superconducting 
side-effects, is observed in (Lai-aCe^) AE. samples with Ce content above x — 
0.0067. However, according to Bader et al. [T^], for x — 0.0064 the expected 
normal-state and measured superconducting-state heat capacities do not differ 
significantly. Thus a mean-field normal-state theory of AC/c for (Lai_ x Ce x ) AE. 
with x — 0.0064 can be reliable. 

The number of valence electrons per host atom in LaAE equals 8/3. For 
ai= 0.0064, 

_ 0.0064 _ _1_ _ 8 _ a _ 3742 

° ~ 2.9936 ~ 1871 ' Ul ~ 3 C ~ 3 ' 

Since the spin of Ce ions equals 1/2 [12] , therefore 

h (0 = irr~j= t anh iPi^i) (32) 
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9 




10 




11 



1,4 




0,0 H 1 1 1 

0,01 0,1 1 10 

T(K) 

Figure 4: The mean-field impurity heat capacity of (Lai-^Ce^) AI2 with x — 
0.0064 and values of 60, &i, 7 5 M, AT given in Table[TJ The points are experi- 
mental results from Ref. [T2"] . 



and AU m /c for a sample S of (Lai^Ce^) AI2, expressed in joules, equals 

c^AL^ = ( i 7 2 - JmV/I (0 - mV/ 2 (£) [6 

^ ^ (33) 
+&1/2® +•••]) I ■ 602 - 2 ' 160 - 2 ' 1871A/" 1 

if Eq. (|44l is used and 60, 7, £ are given in powers of eV. 

The mean- field AC/c curve best fitting to experimental AC/c of Ref. [T2"] 
was obtained for M = 40 and is depicted in Fig. [U The corresponding values 
of other parameters are given in Table [T] Since the error of experimental AC/c 
values is relatively high above 5K [T2], the mean-field AC/c curve in Fig. Q] 
provides a good fit to experiment. 

The equation /s(£) = £ for £(T) has a unique minimising solution with £(T) G 

(ll ■ 10~ 5 VeV, 19 • 10~ 5 VeV) for T E (0.015 K, 10 K). The graph of £(T) is, 
similar as those for CuCr in Fig. [H increasing in T. 
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4 Concluding remarks 



Progress has been made in improving the quality of theoretical impurity heat 
capacity graphs AC(T) of DMA. Scaling procedures, used in previous s-d theo- 
ries to adjust thermodynamic functions to experimental data, were unnecessary 
and nonlinear dependence of AC(T) on impurity concentration c in CuCr has 
been accounted for the first time. 

The mean-field theory of a dilute s-d system established in Ref. [TT] has 
also proved capable of explaining the temperature dependence of impurity heat 
capacity AC of (Lai-zCe^) AI2 with x = 0.0064. Partial quantitive agreement 
between theory and experiment has been achieved for CuCr alloys, especially 
for the smaller value of impurity concentration considered. 

Higher order terms of the expansion of /1 should be included in order to 
improve these results. 



The integral ([15)) is an M -fold product of 3-dimensonal integrals 

r p 3 

/ = J d 3 i?cxp -7^2 J ^2 {R q - r Qq - uj q (s)) 2 ds 
A n 1 =1 



(34) 



and exp [— f3Muo\. The representation (fl"6|) of (p(r, r 0/ A > ror large /?, obtains 

by apllying to the / the dominated convergence theorem (in order to express the 
integral in the exponent as a limit of partial sums) and the method of steepest 
descent: 



I = lim / d 3 i?exp 



^~ u2 £ ( R i - r °i - Xk <^ 

k=l q 



lim / d i?exp 

m— 

.JL 

2m 



2m 



U2 X kq 
kq 



11 2 



lim 



2tt \ 

i^oa \f3ll2 J 



3/2 



exp 



JL 

2m 



kq q \ k J 



\(3u 2 



exp 



P / P 

-\u2 /o; 2 ( s )d S + i/3-VE / 
« \o 



W q (s)ds 



(35) 
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where Xk q = u q y^§J ■ The squared integral in the exponent can be linearized, 
using the identity 



cxp(a 2 ) = (27T)- 1 / 2 J cxp + V2a^ 



di 



which yields 



/ 2tt n 








^ 3Af/2 


exp 


j-s/a ( 


' 2^ X 


3M/2 

) 















/3 / 18 

-^Mu 2 f u: 2 (s)ds + ^- 1 Mu 2 J2 [ 
" 9 Vo 

9 ^ 



w g (s)ds 



d£id£ 2 d£ 3 exp 



^Mb 2 / w 2 (s)ds+ y/p- 1 Mu 2 ^2,^ q J uj q (s)ds 



f 2?r 



|(M-1) 



M 3/2 / dAidA 2 dA 3 exp 



-Mu 2 ^ / (U) q (s) - Xqfds 



(36) 



\/3u 2 

in accord with Eq. (jT6 

B 



The form of /i defined by Eq. J7J) is crucial in the mean-field theory of the 
s-d system. Here fi is derived under the assumption that h^' is de- 

fined in terms of a function —gy/nU (\R a — f]j , which has absolute maximum 
equal «o and absolute minimum equal uq, and that the expansion in (|15p is 
terminated on the first term. The space average of exp 
impurity positions then equals 



0) 



over 



<'xi> -/3^' M) (e,0) 
and Eq. (|10p for 2 takes the form 



dpp 2 < ( z 1 exp 



z exp 



exp 
exp 



forer z = 1 
forer z = — 1 



47TD 



2m 



(3 -Mv {[ + 

2m 



(37) 



14 



where v = |A|n 1 . Following Huang [22], one transforms Eq. (|57|) to 

v\ 3 (/ 3/2 (ze-^«o«) + / 3/2 (ze«)) = 1, (38) 

where 



v\ 3 = (l3e F y 3/2 . 

For large z, fs/2{z) is given by Sommerfeld's expansion: 



f 3/2 (z) = ^= {{\nzf 2 + ^ (lnz)- 1/2 + ...). (39) 

Let Mwo^Ep 1 = /i, Mwo^Ep 1 = v. Then z satisfying Eq. ([35]) is given by the 
expansion 

lnz = aif3£-F + a + a_i (/3£p) _1 + . . . (40) 
with ^ 

»i(0 = 1 + 2 - ") - is (ai + v ? + ' • ■ 
a = 
1 \ -V2 

1 - \2 



«-i(0 = "12 ^-4^ + ") J 12 ' S ( ' 

;From Eq. one hnds (cf. Ref. [TT]) 



Hence 



A(0 = (Afn/3)" 1 ^ Trln (l + 



with 



h(0 = «A 3 (-.0/3/2 (^e-^-S) + « / 3/2 (ze*^)) 
= bo + 6i£ - -3-60 (/3ep) +... 



°o = 7; ( v o - u ) , bi = - — (u + v ) 



(41) 



The energy of the electrons' subsystem also expresses in terms of f±. According 
to Eqs. J3ED, © 
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The first term on the rhs of Eq. (|4"2|) can be evaluated in a similar manner as 
the energy of a free Fermi gas, viz., 

„r"F7-W\ 



^/^{(In^-^tioK-i/)) 1 



+ ^tt 2 (lnz-/3Aftio (e-??)) 1/2 + ... 

+ (In z + /?Mv (£ - ??)) 5/2 + (In z + f3Mv (£ - il)) 1 ' 2 + 

8 

= ^e F 1 2 + 5a_x (£ - rj) (Ps F )- 2 + V (/? £f )" 2 

+ ^M 2 (u + vo fe F 2 (t-v) 2 + - 

= (nT?IT«) 1} - ^nM& a (£ - r,) 2 (/3 £f )- 2 + InMfc (£ - V ) 2 + . . . 

(43) 

If the Taylor expansions of U a (r) around the absolute maximum and absolute 
minimum are cut off on the second derivatives, the structure of /i, as well as 
that of the expectation values ([42]) , ([43]) , changes, e.g. under the approximations 
of Section [21 fi is no longer regular at £ = 0: 

/i(0 = « { Ao (t (« a , /3))- 3 aZo( ^'" 2 ° /3/2 (^o (urf, uaO) 

+ A (t (-V 2 , Pj) 7^ /3/2 (^A) (-«o£, -W 2 0) J . 

with 

v\^(t) = ±3^(te F )- 3/2 

(u2 (fa) denoting the second derivative of —g^/nU (jl a — at the absolute 

minimum (maximum)) and the corresponding coefficient b^, in the new expan- 
sion of /1 around £0 7^ 0, will differ from b'[ resulting in the modified Eq. (j43|) . 
In view of the presumable smallness of 

ACTo = (nr?iT«) -(nr?JJ«) 

where he' M ^ — nT^hi 1 '^, these modifications can be exploited in the simplest 
manner by restricting uq, vq, 112, V2 to values for which b'{ vanishes. As a 
consequence, 

AU a « (44) 

whereas the range of b[ can be expected to include also negative values. Thus, 
in general, 

m) = b' + b^+... , b' Q ,b[ eR 1 (45) 
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